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Abstract 



Q \ The inclusive photoproduction of D* mesons and of L'*-tagged dijets is investigated 

CN ' with the HI detector at the ep collider HERA. The kinematic region covers small photon 

virtualities Q^ < 2 GeV^ and photon-proton centre-of-mass energies of 100 < W-^p < 
285 GeV. Inclusive D* meson differential cross sections are measured for central rapidities 

K> ■ \r]{D*)\ < 1.5 and transverse momenta pt{D*) > 1.8 GeV. The heavy quark production 

5—1 ■ process is further investigated in events with at least two jets with transverse momentum 

- - -' pTiJQt) > 3.5 GeV each, one containing the D* meson. Differential cross sections for 

L'*-tagged dijet production and for correlations between the jets are measured in the range 
\r]{D*)\ < 1.5 and pt{D*) > 2.1 GeV. The results are compared with predictions from 
Monte Carlo simulations and next-to-leading order perturbative QCD calculations. 
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1 Introduction 

At the electron-proton collider HERA charm quarks are predominantly produced via boson 
gluon fusion, 7(7 — ;■ cc, where the photon is emitted from the incoming lepton and the gluon 
originates from the proton. The cross section is largest for photoproduction, i.e. for photons 
with negative four- momentum squared (virtuality) Q^ ~ GeV^. In addition to hard direct 
scattering off the photon, processes have to be considered in which the partonic structure of the 
photon is resolved. The charm quark mass provides a hard scale which justifies the applicability 
of perturbative QCD (pQCD). 

Previous measurements of the photoproduction of charm quarks at HERA cover inclusive 
D* meson production [[13-131, production of D* mesons with associated dijets [[D[3]-[5l and heavy 
quark production using events with a D* meson and a muon [6[. In this paper, single and 
double differential cross sections are presented for the inclusive production of D* mesons and 
the production of two jets with one of the jets containing the D* meson. They are compared 
to leading and next-to-leading order pQCD predictions using different hadronisation models. 
Compared to the previous HI analysis of inclusive D* photoproduction [[3|, a seven times larger 
signal sample is analysed here. 

Studying events in which at least two jets could be reconstructed, with one of the jets con- 
taining the D* meson, allows further investigations of the details of the heavy quark production 
process. The jets are measured down to transverse momenta of prQet) = 3.5 GeV. While 
the jet containing the D* meson originates from a charm or anticharm quark produced in the 
hard subprocess, the non-Z}*-tagged jet, refered to as other jet, can result from either the other 
heavy quark or a light parton (e.g. a gluon). Correlations between the two jets are studied using 
variables which are sensitive to higher order effects and to the longitudinal as well as to the 
transverse momentum components of the partons entering the hard scattering process. 



2 QCD Calculations 

The data presented in this analysis are compared with Monte Carlo simulations based on leading 
order (LO) matrix elements supplemented by parton showers and with next-to-leading order 
(NLO) calculations. The calculations are performed using either the coUinear factorisation or 
the A;t-factorisation approach. The coUinear factorisation makes use of the DGLAP [7] evolution 
equations, while in fct-factorisation the CCFM [8] evolution equations are employed. In the 
coUinear approach transverse momenta obtained through the initial state QCD evolution are 
neglected and the transverse momenta are generated in the hard scattering process. Effects 
from the non- vanishing transverse momenta of the gluons enter only at the NLO level. In the 
fcf-factorisation ansatz the transverse momenta of incoming gluons, /c^, are already included 
at leading order both in the off-shell matrix element and the /cj-dependent unintegrated gluon 
density [|9|. Corrections appearing only at higher order in coUinear factorisation are hence 
partially included at LO in the fcj-factorisation approach. 

For charm quark photoproduction two classes of processes occur, the direct-photon and the 
resolved-photon processes. In the direct processes the photon emitted from the beam lepton 



enters directly the hard interaction, whereas in the resolved processes the photon acts as the 
source of incoming partons, one of which takes part in the hard interaction. The distinction 
between these two classes depends on the factorisation scheme and the order in which the 
calculation is performed. 

The production of heavy quarks is calculated either in the massive scheme, where heavy 
quarks are produced only perturbatively via boson gluon fusion, or in the massless scheme, 
where heavy quarks are treated as massless partons. These two schemes are expected to be 
appropriate in different regions of phase space [10]: the massive scheme is expected to be 
reliable when the transverse momentum px of the heavy quark is of similar size compared 
to the charm mass nic, whereas the massless scheme is expected to be valid for px ^ rric. 
In the general-mass variable-flavour-number scheme (GMVFNS) a smooth transition from the 
massive to the massless scheme is provided. The structure of the proton and of the photon are 
described by parton distribution functions (PDFs), that have been determined by fits to data in 
various heavy flavour schemes and at different orders of pQCD . 

Monte Carlo (MC) generators are used to simulate detector effects in order to determine the 
acceptance and the efficiency for selecting events and to estimate the systematic uncertainties 
associated with the measurement. The generated events are passed through a detailed simulation 
of the detector response based on the GEANT simulation programm [llj and are processed 
using the same reconstruction and analysis chain as is used for the data. The following two MC 
generators are used: 

PythiA: The MC program Pythia [[13 is based on LO QCD matrix elements with leading- 
log parton showers in the coUinear factorisation approach. Pythia includes both direct 
photon gluon fusion and resolved-photon processes. In the resolved-photon processes 
either a charm quark or a gluon from the photon enters the hard scattering. In the inclusive 
mode of Pythia used here charm quarks are treated as massless partons in all steps of the 
calculation in both types of processes. The hadronisation process is simulated using the 
Lund string fragmentation model [[T3l . The Bowler fragmentation model |[T4l is applied 
to fragment the charm quark into a D* meson. The longitudinal part of the fragmentation 
is reweighted to the parameterisation by Kartvelishvili et al. [l5j which depends on a 
single parameter a. The latter is set to the values determined by HI [16J, which depend 
on the centre-of-mass energy squared of the hard subprocess s (see table [I])- The proton 
structure is described by the PDF set CTEQ6L [17J. For the photon the PDF set GRV-G 
LO |[18J is used. 



Cascade: The Cascade [[T9|| MC program is used for simulating events based on LO QCD 
calculations in the /cj-factorisation approach. The direct boson gluon fusion process is 
implemented using off-shell matrix elements and incoming gluons which can have non- 
vanishing transverse momenta. Higher order QCD corrections are simulated with initial 
state parton showers applying the CCFM evolution [8]. The unintegrated PDFs of the 
proton from set AO [|2P1 are used. The hadronisation of partons is performed with the 
Lund string model as implemented in Pythia. For the fragmentation of the charm quarks 
into D* mesons the same reweighting procedure to the parameterisation of Kartvelishvili 
et al. is applied as in the case of Pythia. 



Fragmentation parameter a. 






Pythia 


Cascade 




^threshold 

[GeV^] 


a for 

■5 < Sthreshold 


a for 

■5 > Sthreshold 


a for 

S < Sthreshold 


a for 

S > Sthreshold 


Central value 


70 


10.3 


4.4 


8.4 


4.5 


Variations 


70 
70 
50 
90 


8.7 
12.2 
10.3 
10.3 


3.9 

5.0 

4.4 
4.4 


7.3 

9.8 

8.4 
8.4 


3.9 

5.1 
4.5 
4.5 



Table 1: Fragmentation parameters a in the Kartvelishvili parameterisation used in the MC 
simulations. In the two regions of the invariant mass squared of the cc pair, s, separated by the 
boundary Sthreshold, two different values of a are used. 



For the comparison of data with NLO predictions, calculations based on the massive ap- 
proach and the general mass variable flavor number scheme are used. The uncertainties of the 
calculations are estimated by varying the charm mass, rric, the factorisation scale, /i/, and the 
renormalisation scale, fir- The detailed settings are given in tabled For the comparison in the 
D*-tagged dijet sample only MC@NLO is used since it provides a full hadronisation of the 
final state. 





FMNR 


GMVFNS 


MC@NLO 


Parameter 


Central 


Variations 


Central 


Variations 


Central 


Variations 


Charm mass nic/GeV 
Renorm. Scale fir/mx 
Fact. Scale fif/mx 


1.5 
1 
2 


1.3 1.7 

0.5 2 

1 4 


1.5 
1 
1 


0.5 2 
0.5 2 


1.5 
1 
2 


1.3 1.7 
0.5 2 

1 4 



Table 2: Parameters and variations used in the NLO calculations of FMNR 1127112211 . 
GMVFNS t26.^ and MC@NLO ifSTIl . 



FMNR: The FMNR program [[211122] is based on an NLO calculation in the massive scheme 
in the coUinear approach. The resolved and direct processes are calculated separately. 
The program provides weighted parton level events with two or three outgoing partons, 
i.e. a charm quark pair and possibly one additional light parton. The fragmentation of a 
charm quark to a D* meson is treated by a downscaling of the three-momentum of the 
quark in the charm-anticharm rest frame according to the Peterson fragmentation function 
with a parameter value of e = 0.035. The PDF sets HERAPDFl.(U[[23l for the proton 



'The HERAPDFl.O set was determined from inclusive deep-inelastic scattering data from the HI and ZEUS 
experiments in the GMVFNS. It has been checked that the difference to a PDF set determined in the massive 
scheme, CTEQ5F3 |24|, is significantly smaller than the effect of the variations considered for the systematic 
uncertainty of the FMNR predictions. 



and GRV-G HO IfTSlI for the photon are used. For the strong coupling, the five-flavour 
QCD scale Ag^^ is set to 0.2626 GeV. The charm mass is set to rric = 1.5 GeV and 
varied by ±0.2 GeV for an uncertainty estimate. This variation covers the central value 
for the pole mass of the charm quark [25 J. The renormalisation and factorisation scale 
are set to fi,- = rriT and fif = 2 ■ rriT with mx being the transverse massdefined as 
m^ = ml + {p\ ^ + p\^/2, with p^^c and px^c denoting the transverse momenta of the 
charm and anticharm quark, respectively. In order to estimate the uncertainties related to 
missing higher orders, the renormalisation and factorisation scales are varied by a factor 
2 up and down. Each variation is done independently, leading to in total 6 variations. 
The resulting uncertainties are added in quadrature separately for positive and negative 
deviations to obtain the total uncertainties. 

GMVFNS: A next to leading order cross section prediction for direct and resolved contribu- 
tions to the cross section has been provided in the GMVFNS [I26U27 1. The transition from 
the charm quark to the D* meson is given by the KKKS fragmentation function which 
takes DGLAP evolution and finite-mass effects into account |[28l . The parton contents 
of the proton and of the photon are described by the PDF sets HERAPDFl.O [23J and 
AFG04 [29|, respectively. The charm mass is set to rric = 1.5 GeV, and the renormalisa- 
tion and factorisation scales are chosen to be n,- = fif = rriT- The uncertainties related to 
missing higher orders are estimated by varying the renormalisation scale, the factorisation 
scale for the initial state and the factorisation scale for the final state independently by a 
factor 2 up and down while satisfying the condition that the ratio of any of the two scales 
is 1/2, 1 or 2. This leads to 14 independent variations. The maximum and minimum 
values found by this procedure are used to determine the systematic uncertainty [|27l . 



MC@NLO: In the MC@NLO framework fSOl, predictions for heavy flavour production at 
HERA [31 J are provided which combine an NLO calculation in the massive approach 
with parton showers and hadronisation. The direct and resolved part of the cross sec- 
tion are calculated separately. MC@NLO uses parton showers with angular ordering 
to simulate higher order contributions and the cluster fragmentation as implemented in 
HERWIG |32|. A factor of 1.34 is applied to the MC@NLO predictions in order to cor- 
rect the c — )■ D* branching fraction in HERWIG to the experimental value [33 J. The PDF 
sets HERAPDFl.O [23J for the proton and GRV-G HO [18J for the photon are used. For 
an estimation of the uncertainty, the charm mass and the renomalisation and factorisation 
scales are varied separately, and the resulting uncertainties are added in quadrature. 



3 HI Detector 



A detailed description of the HI detector can be found elsewhere [|34l . Only the components 
essential to the present analysis are described here. 

The origin of the HI coordinate system is the nominal ep interaction point. The positive 
z-axis (forward direction) is defined by the direction of the proton beam. Transverse momenta 
are measured in the x-y plane. Polar {6) and azimuthal {ip) angles are measured with respect to 
this reference system. The pseudorapidity is defined as r^ = — lntan(6'/2). 



Charged particles are measured within the central tracking detector (CTD) in the pseudo- 
rapidity range —1.1 A < r] < 1.1 A. The CTD comprises two large cylindrical jet chambers 
(inner CJCl and outer CJC2) and the silicon vertex detector [[351 . The CJCs are separated 
by a drift chamber which improves the z coordinate reconstruction. A multiwire proportional 
chamber mainly used for triggering [36] is situated inside the CJCl. These detectors are ar- 
ranged concentrically around the interaction region in a solenoidal magnetic field of 1.16 T. 
The trajectories of the charged particles are measured with a transverse momentum resolution 
of a{pT)/pT ~ 0.5%pT/GeV © 1.5% [i37il . The CJCs also provide a measurement of the 
specific ionisation energy loss dE/dx of charged particles. The interaction vertex is recon- 
structed from CTD tracks. The CTD also provides trigger information based on track segments 
measured in the CJCs [|38l . At the first two levels of this fast track trigger (FTT) tracks are re- 
constructed online from the track segments in the CJCs. At the third level of the FTT invariant 
masses of combinations of tracks are calculated [i39ll40l . 

Charged and neutral particles are measured with the liquid argon (LAr) calorimeter, which 
surrounds the tracking chambers. It covers the range —1.5 < r] < 3.4 with full azimuthal 
acceptance. Electromagnetic shower energies are measured with a precision of a{E)/E = 
12%/^/EjGeV © 1% and hadronic energies with a{E)/E = 50%/y/E/GeV © 2%, as de- 
termined in test beam measurements [HD . A lead-scintillating fibre calorimeter (SpaCal) [42] 
covering the backward region —4.0 < i] < —1.4 completes the measurement of charged and 
neutral particles. For electrons a relative energy resolution of a{E)/E = 1%/ ^jE/GeY © 1% 
is reached, as determined in test beam measurements ^^i\. The hadronic final state is recon- 
structed using an energy flow algorithm which combines charged particles measured in the CTD 
with information from the SpaCal and LAr calorimeters [04]. 

The luminosity determination is based on the measurement of the Bethe-Heitler process 
ep — )■ ej97 where the photon is detected in a calorimeter located dX. z = —104 m downstream of 
the interaction region in the electron beam direction. 



4 Event Selection and Reconstruction 

The data sample was recorded in the years 2006 and 2007, when electrons with an energy of 
27.6 GeV were collided with protons with 920 GeV. 

Photoproduction events are selected by requiring that no isolated high energy electromag- 
netic cluster, consistent with a signal from a scattered electron, is detected in the calorimeters. 
This limits the photon virtuality to Q^ < 2 GeV^. 



4.1 Inclusive D* Sample 

The triggering of the events relies on the reconstruction of the final state particles originating 
from the D* decay. For this purpose all three levels of the FTT are used. At the first level, 
where tracks are reconstructed only in the transverse plane, the selection criteria are based 
on track multiplicities above certain transverse momentum thresholds. These conditions are 
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refined on the second level, and on the third level invariant masses and charge combinations 
consistent with the decay channel D*^ — )■ D'^ti^^^^ — )• K^Ti^iifi^^ are required [|40| . Three 
trigger conditions with different thresholds for the transverse momentum of the D* candidate 
are used. The analysis is therefore performed in three separate pt{D*) regions corresponding 
to the different luminosities: C = 30.7 pb"^ for 1.8 < priD*) < 2.5 GeV, C = 68.2 pb"^ 
for 2.5 < PriD*) < 4.5 GeV, and C = 93.4 pb"^ for priD*) > 4.5 GeV. The requirement 
that all decay particles have to be in the acceptance of the CJC limits the analysis to central 
rapidities for the D* meson \'r]{D*)\ < 1.5 and photon-proton centre-of-mass energies in the 
range 100 < W^p < 285 GeV. 



The 7p centre-of-mass energy is reconstructed using the Jacquet-Blondel method P3I : 
W^7P = y/yjB s with yjB = ^HFsi-^ ~ Pz)i/{'^ Ee), whcrc s and E^ denote the square of 
the ep centre-of-mass energy and the energy of the incoming electron, respectively, and the sum 
'^HPS runs over the energy E and the longitudinal momentum pz of all final state particles. 
The D* inelasticity z{D*), which corresponds to the fraction of photon energy transferred to 
the D* meson in the proton rest frame, is defined by z{D*) = P ■ p{D*)/{P ■ q), with P, p{D*) 
and q denoting the four- momenta of the incoming proton, the D* meson and the exchanged 
photon, respectively. It is reconstructed as z(D*) = (E — Pz)d*/{'^ Vjb E^,). The inelastic- 
ity distribution is sensitive to the kinematics of the production mechanism and to the c —)■ D* 
fragmentation function. 

The D* meson is detected via the decay channel D*^ — )■ D'^nfi^^ — )■ K^TT^T^fioy^ with a 
branching fraction of BTZ = 2.63 ± 0.04% [|25l . The tracks of the decay particles are recon- 
structed using the CTD information. The invariant mass of the K^tt"^ system is required to be 
consistent with the nominal D^ mass [25J within ±80 MeV. The signal to background ratio is 
improved by applying a loose particle identification criterion to the kaon candidates based on 
the measurement of the specific energy loss, dE/dx, in the CTD. In addition the background is 
reduced by a cut on the fraction of the transverse momentum carried by the D* with respect to 
the scalar sum of transverse energies of the hadronic final state, excluding the forward region 
(6 < 10°). This fraction is required to be PT{D*)/{^j^pg Et^i) > 0.1. This criterion accounts 
for the harder fragmentation of charm compared to light flavours. 

The D*^ candidates are selected using the mass difference method ['461. In figure [T^) the 
distribution of the mass difference AM = m^Kmisiow) — rn^Kix) of the final D* candidates is 
shown. A clear peak is observed around the nominal value of AM = 145.4 MeV dlSl. 

The wrong charge combinations, defined as K^-k^-k^i^ with K^-k^ pairs in the accepted 
D^ mass range, are used to constrain the shape of the combinatorial background in the signal 
region. The number of reconstructed D* mesons N(D*) is extracted in each analysis bin by a 
log-likelihood fit simultaneously to the right charge and the wrong charge AM distribution. For 
the signal which has a tail towards larger AM values the asymmetric Crystal Ball function [|47l 
is used. The shape of the background is parametrised with the Granet function [i48il . The fit 
is performed in the RooFit framework ^49\. The fit to the inclusive data sample yields 8232 ± 
164 D* mesons. To improve the convergence of the fit in each analysis bin, the parameters 
describing the asymmetry of the Crystal Ball function are fixed to the values found by the fit to 
the complete data set. The width of the peak varies in dependence on the D* kinematics and is 
therefore left free. More details can be found in [JSOl. 



4.2 D* -tagged dijet Sample 

For the selection of the D* meson in the D*-tagged dijet sample, the requirements are the same 
as for the inclusive D* sample, except that the requirement on the specific energy loss dE/dx is 
removed, and the cut on pt{D*) is increased to 2.1 GeV because of large backgrounds at small 
transverse momenta. 

Jets are defined by the inclusive A;i-algorithm [51 J in the energy recombination scheme with 
jet size AR = ^/(Ar])^ + {Aipy = 1 where Aip is expressed in radians. The jet algorithm 
is applied in the laboratory frame to all reconstructed particles of the hadronic final state. To 
prevent the decay particles of the D* candidate from being attributed to different jets, the D* 
candidate is used as a single particle in the jet algorithm, replacing its decay products. In this 
way the jet containing the D* meson {D* jet) is unambiguously defined for each D* candidate. 
In events which contain more than one D* candidate, the jet algorithm is run separately for 
each candidate, and all candidates for which the dijet selection criteria are fulfilled enter the 
AM distribution. The pseudorapidity of the D* jet is restricted to the same range as is used 
for the D* meson, \r]{D* jet)| < 1.5. In addition to the D* jet a second jet is required. Both 
jets have to satisfy prOet) > 3.5 GeV. If there is more than one jet that does not contain the 
D* meson, the one with the highest prOet) is chosen as the other jet. The pseudorapidity of 
the other jet has to be in the range —1.5 < r7(otherjet) < 2.9. The invariant mass Mjj of 
the D* jet and the other jetis required to satisfy Mjj > 6 GeV in order to select jets from the 
partons originating from the hard interaction. More details on the selection of the D*-tagged 
dijet sample can be found in [52,|. 

The number of -D*-tagged dijet s is extracted from the AM distribution of the D* candidates 
with the same procedure as used for the inclusive D* measurement. The AM distribution for 
the selected events in the dijet sample is shown in figure [Ih). The fit yields asignalof 3937±114 
D* mesons. 

The kinematic range of the inclusive D* measurement and of the D*-tagged dijet measure- 
ment are summarised in table |3] 



5 Cross Section Determination and Systematic Errors 

The bin averaged visible differential cross section with respect to a variable Y (with bin width 
AY) is calculated according to 

da^isjep ^ e D* + X) _ N{D*){l-r) 

dY AY -C-BK-e ^ ^ 

where C is the integrated luminosity, BIZ is the branching ratio of the analysed decay chain 
£)*± _v. jjo^±^^ _^ j^=fj^±^±^^ and (1 — r) a correction factor to account for reflections from 
other D^ decays. The efficiency e includes the detector acceptance, trigger and reconstruction 
efficiencies and migrations between bins. The contributions of D* mesons originating from 
beauty production and from gluon splitting from light flavour production is not subtracted. It is 
estimated from MC predictions to be below 2%. 
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inclusive D* meson and D* -tagged dijet production 


Photon virtuality 


g2 < 2 GeV^ 


7P centre-of-mass energy 


100 < W^p < 285 GeV 


Pseudorapidity of D*^ 


\r]{D*) < 1.5 


inclusive D* meson production 


Transverse momentum of D*^ 


Pt{D*) > 1.8 GeV 


D* -tagged dijet production 


Transverse momentum of D*^ 


PriD*) > 2.1 GeV 


Transverse momentum of D* jet 


PT{D*iet) > 3.5 GeV 


Pseudorapidity of D* jet 


\v{D*ict) <1.5 


Transverse momentum of other jet 


PT(other jet) > 3.5 GeV 


Pseudorapidity of other jet 


-1.5 < r7(other jet) < 2.9 


Dijet invariant mass Mjj 


Mjj > 6 GeV 



Table 3: Definition of the kinematic range of the measurements. 

The systematic uncertainties are determined in each bin separately and are summarised in 
tablelUfor the total cross section. They are divided into uncertainties which are considered to be 
uncorrelated between the bins and uncertainties which change the cross section normalisation 
in all bins. The numbers for the uncertainties listed below are given in per cent of the cross 
section values. 

The following uncorrelated systematic uncertainty sources are considered: 

Trigger Efficiency: The simulation of the FTT is verified by a comparison to data in a sample 
of D* mesons in deep-inelastic scattering triggered by the scattered electron. For the total 
inclusive D* sample the efficiency agrees within a relative uncertainty of 7.5%. This is 
one of the dominant systematic uncertainties. For the D*-tagged dijet sample the trigger 
efficiency is higher, leading to a smaller uncertainty of 3.1% for the total cross section. 

Signal Extraction: For the determination of the uncertainty of the signal fit, different param- 
eterisations for the signal and background functions are used. The resulting uncertainty 
amounts to 1.5%. 

D^ mass cut: The loss of D* mesons due to the D° mass cut is compared between data and 
simulation as a function of the D* transverse momentum, assuming a Gaussian resolution 
for the D^ mass reconstruction. They agree within 2%, which is assigned as uncertainty. 

Reflections: The amount of reflections r from decay modes of the D^ meson other than Z5" — > 
K^tt"^ amounts to 3.8% in the simulation [53| . It is independent of kinematic quantities 
within 1%, which is used as systematic uncertainty. 

Background from deep inelastic scattering: The background originating from deep inelastic 
scattering events is estimated with the RAPGAP [54] MC generator. It is found to be 
below 1%, which is not subtracted but treated as an uncertainty. 
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dE/dx cut: The efficiency of the cut on the dE/dx likelihood of the kaon candidate is studied 
for data and MC simulation in bins of the transverse momentum of the D* meson. The 
relative difference of 1.5% is corrected for in the MC sample. An uncertainty of 0.5% is 
assigned, covering the possible pt{D*) dependence of this correction. 

Hadronic energy scale: The energy scale of the hadronic final state has an uncertainty of 2% 
leading to an uncertainty of the cross section of 0.6% in the inclusive D* sample and of 
2.0% in the D*-tagged dijet sample. 

Model: For the determination of the cross section the Pythia and CASCADE predictions are 
reweighted to describe the data distributions where necessary. For the correction of the 
data the efficiency from the Pythia MC is used. The difference to the efficiency from 
Cascade is taken as a systematic uncertainty. It amounts to 2% (1.5%) for the total 
inclusive D* (Z}*-tagged dijet) cross section. 

Fragmentation: The a parameter of the Kartvelishvili function and the position of the s thresh- 
old are varied within the values given in table [Hresulting in an uncertainty of 2.5% (2.0%) 
for the total inclusive D* (D*-tagged dijet) cross section. 

The following normalisation uncertainties are considered: 

Track finding efficiency: The systematic uncertainty on the track efficiency of 4. 1% per D* me- 
son arises from two contributions: (i) The comparison of the track finding efficiency in 
data and simulation leads to an uncertainty of 2% for the slow pion track and 1% for the 
tracks of the D'^ decay particles, and the uncertainty is assumed to be correlated between 
the decay particles; (ii) the efficiency with which a track can be fitted to the event vertex 
leads to a systematic error of 1% per D* meson. The uncertainty on the track finding 
efficiency is considered to be half correlated between the bins of the measurement. 

Luminosity: The uncertainty on the luminosity measurement for the data sample used in this 
analysis amounts to 5%. 

Branching Ratio: The uncertainty due to the D* branching ratio is 1.5% ||25l . 

All sources of systematic errors are added in quadrature resulting in a systematic uncertainty 
of 10.9% (8.5%) for the total cross section of the inclusive D* (D*-tagged dijet) production. 



6 Results for Inclusive D* Meson Production 

The total visible cross section for D* meson photoproduction is measured to be: 

(^visiep ^eD* + X) = 41.1 ± 0.8 (stat.) ± 3.6 (unc.sys.) ± 2.7 (norm.) nb (2) 

in the kinematic range defined in table |3] The corresponding predictions from Pythia and 
Cascade amount to 43.7 nb and 32.9 nb, respectively. Due to the fact that these predictions 
are based on leading order matrix elements the uncertainty on the normalisation of the cross 
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Uncertainty source 


D* 


D*-tagged dijet 


Uncorrelated uncertainties 


Trigger efficiency 


7.5% 


3.1% 


Signal extraction 


1.5% 


1.5% 


D^ meson mass cut 


2.0% 


2.0% 


Reflections 


1.0% 


1.0% 


Background from deep-inelastic scattering 


1.0% 


1.0% 


dE/dx cut 


0.5% 


— 


Hadronic energy scale 


0.6% 


2.0% 


Model 


2.0% 


1.5% 


Fragmentation 


2.5% 


2.0% 


Track finding efficiency (half) 


2.9% 


2.9% 


Total uncorrelated 


9.2% 


6.0% 


Normalisation uncertainties 


Track finding efficiency (half) 


2.9% 


2.9% 


Luminosity 


5.0% 


5.0% 


Branching ratio 


1.5% 


1.5% 


Total normalisation 


6.0% 


6.0% 


Total 


10.9% 


8.5% 



Table 4: Summary of all sources of systematic uncertainties and their effect on the total D* 
and the £)*-tagged dijet production cross section with the breakdown into sources leading to 
bin-to-bin uncorrelated uncertainties and sources leading to normalisation uncertainties. 



sections is large, and is not quantified here. The NLO calculations predict 26 t^l nb for FMNR, 
37 til nb for GMVFNS and 30 1? for MC@NLO. 

The measured single differential cross section as a function of the transverse momentum 
Pt{D*) and the pseudorapidity r]{D*) of the D* meson, the photon-proton centre-of-mass 
energy Wjp and D* inelasticity z{D*) are presented in table |5]and in figures |2]and |3] The 
data are compared to Pythia, Cascade and the NLO predictions of FMNR, GMVFNS and 
MC@NLO. Since all the predictions have large normalisation uncertainties, the normalised ra- 
tio R^°^^ of theory to data is shown in order to compare the shape of the various predictions to 

the data. R^°^^ is defined as 

1 da^''^^ 



K" 



a. 



calc 



dy 



da' 



data 



(3) 



a[ 



data 



dY 



where a^\l^ {cr'^f^'^) and da'^^'^/dF (da'^^^^/dF) are the total and differential cross section of the 
model under consideration (of the data), respectively, and Y denotes any measured variable. 
In this ratio the normalization uncertainties of the data (luminosity, branching ratio and half of 
the tracking uncertainty) cancel. Similarly, uncertainty sources of the NLO predictions altering 
the normalisation only do not affect R^°^^ since for each variation the total and the differential 
cross section are varied simultanously. 
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The single differential cross sections are compared to the predictions of the LO MC simula- 
tions in figure |2l The steep decrease of the cross section with increasing transverse momentum 
Pt{D*) is reasonably reproduced by Pythia, while CASCADE falls slightly slower than the 
data. Both MC simulations describe the shape of the observed r](D*) distribution within un- 
certainties. The cross section decreases as a function of the 7p centre-of-mass energy W-yp, as 
expected from the photon flux in the equivalent photon approximation ||55l . CASCADE predicts 
a smaller fraction of D* mesons being produced at small inelasticities z{D*), similar to what 
has been observed in deep inelastic scattering at HERA [53 1. All distributions are reasonably 
well described by Pythia. 

A comparison of the single differential cross sections to the predictions of the NLO cal- 
culations is shown in figure |3] For all measured quantities the precision of the measurement 
presented here is much better than the estimated uncertainty of the NLO calculations. The un- 
certainty of the NLO predictions is dominated by the variation of the renormalisation scale fir, 
which has a large effect on the absolute cross section, while the differences in the shapes tend 
to be smaller. Within these large theoretical uncertainties, both the FMNR and GMVFNS pre- 
dictions agree with the measured cross section as a function of pt{D*), while the MC@NLO 
underestimates the data at small pt{D*). The pt{D*) shape is best described by the GMVFNS 
calculation, while FMNR and MC@NLO predict a harder spectrum than observed in data as 
can be seen in the ratio i?"""^™. The underestimation of the low pt{D*) region by the central 
FMNR and MC@NLO predictions results in a low normalisation in the other distributions. The 
shape of the i]{D*) distribution is reasonably well described by all NLO calculations. All three 
NLO calculations give a rather precise prediction of the shape of the W-yp distribution, which 
describes the measurement. Given the large uncertainties the predictions for the z{D*) distribu- 
tion agree with the data, although when using the central parameter settings for the calculations 
they differ in shape with respect to data. 

Previous HI and ZEUS analyses of D* meson photoproduction [[Tl|3l, albeit in different 
kinematic ranges in the photon virtuality Q^ and the photon-proton centre-of-mass energy W.yp, 
lead to similar conclusions: while all predictions give a good description of the W^p distribution, 
differences between data and theoretical predictions are observed for variables sensitive to the 
quantities of the outgoing charm quark. 

In order to investigate the correlation between pseudorapidity and transverse momentum, a 
double differential measurement in pt{D*) and r7(D*) is performed (table|6]). The cross sections 
of the leading order MCs Pythia and Cascade in the three Pt{D*) regions shown in figure|4] 
reflect the different pt{D*) dependences seen in figured Both models are in broad agreement 
with the data. The comparison of the NLO calculations with the data in figure |5]leads to similar 
conclusions as for the LO MC programs. 



7 Results for D* Tagged Dijet Production 

The integrated D*-tagged dijet cross section in the visible range given in table[3]is measured to 
be 

o-«s(ep-^ e D* jet + other jet + X) =9.68±0.28 (stat.)±0.51 (unc.sys.) ±0.64 (norm.) nb. 

(4) 
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The corresponding predictions from Pythia, Cascade and MC@NLO amount to 8.9 nb, 
8.1 nb and 7.1 'tfl nb, respectively. In the common range of transverse momentum, pt{D*) > 
2.1 GeV, the ratio of the D*-tagged dijet to the inclusive D* cross section is 0.304 ± 0.013 ± 
0.031, compared to 0.271 and 0.311 for Pythia and Cascade, respectively. MC@NLO pre- 
dicts a ratio of 0.309 lJ][J^i^. 

The bin averaged differential cross section for the D*-tagged dijet production as a function 
of the transverse momentum px and the pseudorapidity 7] of both the -D* jet and the other jet are 
listed in table|7]and shown in figures |6] and |7] On average, the other jet is more forward than the 
D* jet not only due to the larger measurement range in r], but also within the common region of 
—1.5 < T] < 1.5. This behaviour is consistent with the expectation that the other jet originates 
not always from a charm quark. This observation confirms the result of the previous HI analysis 
of D*-tagged dijet photoproduction [[31 with improved precision. In figure [6]the measurements 
are compared to the Pythia and the Cascade predictions. The shapes of the distributions are 
described well by both models. In figure |7] the measurements are compared to the predictions 
of MC@NLO. At low transverse momenta of both the D* jet and the other jet, the predictions 
lie significantly below the measurement. This results in a smaller total visible cross section 
which is also observed in the t] distribution. The uncertainty band of the MC@NLO prediction 
includes both variation of the charm mass and variations of the factorisation and renormalisation 
scales as described in section |2l 

In order to investigate further the charm production dynamics, several variables related to 
the structure of the hadronic final state are studied. The correlation between the jets in the longi- 
tudinal and transverse directions is experimentally assessed by the difference in pseudorapidity 
Ai] = //(other jet) — 7](D* jet) and in the azimuthal angle \A.ip\ between the D* jet and the 
other jet. The amount of QCD radiation in addition to the the two leading jets is investigated 
using the mass variable Mx = a/(-P + q — (ji + J2)Y with P, q, ji and J2 being the four- 
vectors of the initial proton, the exchanged photon, the D* jet and the other jet, respectively. 
In direct photon processes without radiation, Mx is expected to be close to the proton mass, 
whereas resolved processes as well as additional QCD radiation will increase Mx- The fraction 
x^ of the longitudinal photon momentum entering the hard scattering process can be used to 
distinguish direct and resolved processes: in coUinear factorisation at LO a resolved photon 
process is characterised by x^ < 1, while a direct process has x^ = 1. In the Z}*-tagged dijet 
sample, x^ is approximated by 



EjetsJE-P. 



^ ^^et.v ^^.^ ^3^ 



The sum in the numerator runs over the particles in the two selected jets, whereas the sum in 
the denominator contains all reconstructed particles of the hadronic final state. 

In table [8] and figures [8] and [9] the bin averaged differential cross sections for the D* -tagged 
dijet production as a function of the difference in pseudorapidity Ar/ and in azimuthal angle 
|A(/j| between the other jet and the D* jet, the mass Mx and x-y are presented. The cross 
section as a function of Ar] is not symmetric because the other jet is on average more forward 
than the D* jet. The shape in Ar] is reasonably well described by all QCD calculations. The 
cross section as a function of \Alp\ shows a significant contribution away from the back-to- 
back configuration at \Aip\ ~ 180°. Such a configuration can be described by models which 
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include significant contributions from higher order QCD radiation or a transverse momentum 
of the gluon in the initial state. Whereas Pythia predicts a too small relative contribution of 
these configurations, CASCADE overestimates them. The prediction from MC@NLO, shown 
in figure I9J5), agrees well in shape with the measurement. 

The cross section as a function of the invariant mass Mx is reasonably well described by the 
predictions of CASCADE and Pythia in the region of Mx < 120 GeV, whereas the measured 
cross section is larger than the predictions for the highest Mx bin. The large Mx region is 
correlated with the region of small Xj, where also the predictions are below the measurement. 
MC@NLO predicts a different shape for Mx and is not able to describe the shape of the x-y 
distribution. 

The I Ac/9| dependence of the cross sections in two regions of x-y is presented in tableland 
in figure [lOl Pythia is in agreement with the data. Cascade overestimates the contribution 
from small |Ac/9| in both Xj regions. MC@NLO describes the shape well in the region of 
small x^, where resolved photon processes are enhanced, but is too low in normalisation. At 
large x^ values MC@NLO predicts the size of the cross section correctly, but overestimates the 
contribution from small \Aip\. 

The cross sections for D*-tagged dijet production show that in general both hard partons in 
the final state can be described reasonably well by the QCD predictions, while the details and 
especially the correlations between the D* jet and the other jet are not described very well by 
these theoretical calculations. 



8 Conclusions 



The production of Z^* mesons in the photoproduction regime is investigated with the HI detec- 
tor at HERA with a seven times larger signal sample compared to the previous HI measurement. 
The events containing D* mesons were triggered by the tracks of the decay particles in the chan- 
nel D^ — )■ D^T^fiou, ~^ ^^'^'^'^fiow- Single and double differential cross sections are measured, 
and the results are compared to leading order QCD models provided by the MC simulation 
programs Pythia and Cascade and to the next-to-leading order pQCD calculations FMNR, 
GMVFNS and MC@NLO. The precision of the cross section measurements far exceeds the 
predictive power of the NLO theories. The shapes of the differential cross sections, however, 
are less sensitive to the theoretical uncertainties, and generally show reasonable agreement with 
the data. 

The cross section for D*-tagged dijet production is measured and compared to predictions 
of Pythia, Cascade and MC@NLO. The resuks are consistent with the expectation that the 
non-Z}*-jet can originate not only from a charm quark but also from a light parton. Significant 
contributions from higher order QCD radiation or transverse momenta of the partons in the 
initial state are needed to describe the cross section away from the back-to-back configuration 
between the D* jet and other jet at |Ay9| ^ 180°. The cross sections as a function of the 
transverse momentum and the pseudorapidity of the D* jet and the other jet are reasonably well 
described by the predictions. However, significant differences are observed in the description 
of some variables related to the structure of the hadronic final state, such as | Ayj], Mx and x^. 
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HI inclusive D* cross sections 


Pt{D*) range 


d(r/dpT{D*) 


Stat. 


sys. 


[GeV] 


[nb/GeV] 


[%] 


[%] 


1.8 2.1 


36 


±12 


±13 


2.1 2.5 


29 


±8 


±13 


2.5 3.0 


15 


±5 


±11 


3.0 3.5 


8.6 


±6 


±8 


3.5 4.5 


4.3 


±3 


±8 


4.5 5.5 


2.3 


±4 


±9 


5.5 6.5 


0.89 


±5 


±7 


6.5 9.0 


0.25 


±6 


±8 


9.0 12.5 


0.047 


±12 


±11 


ri{D*) range 


da/dri{D*) 


Stat. 


sys. 




[nb] 


[%] 


[%] 


-1.5 -1.0 


13 


±5 


±10 


-1.0 -0.5 


16 


±4 


±10 


-0.5 0.0 


18 


±4 


±10 


0.0 0.5 


15 


±4 


±10 


0.5 1.0 


12 


±5 


±10 


1.0 1.5 


7.9 


±10 


±10 


W^p range 


da/d{Wyp) 


Stat. 


sys. 


[GeV] 


[nb/GeV] 


[%] 


[%] 


100 140 


0.34 


±3 


±10 


140 180 


0.29 


±3 


±10 


180 230 


0.19 


±4 


±10 


230 285 


0.11 


±6 


±10 


z{D*) range 


da/d{z{D*)) 


Stat. 


sys. 




[nb] 


[%] 


[%] 


0.00 0.10 


45 


±14 


±11 


0.10 0.20 


89 


±5 


±11 


0.20 0.35 


76 


±3 


±10 


0.35 0.55 


55 


±3 


±9 


0.55 1.00 


13 


±4 


±11 



Table 5: Bin averaged single differential cross sections for inclusive D* production in bins of 
Pt{D*), r](D*), Wjp and z{D*) with their statistical and uncorrelated systematic uncertainties. 
The normalisation uncertainty of 6.0% is not included. 
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HI inclusive D* cross sections 


1.8 <pt{D*) <2.5GeV 


r]{D*) range 


d'^a/drjdpT stat. 


sys. 




[nb/GeV] [%] 


[%] 


-1.5 -1.0 


13 ±12 


±14 


-1.0 -0.5 


12 ±12 


±14 


-0.5 0.0 


14 ±11 


±13 


0.0 0.5 


10 ±16 


±13 


0.5 1.5 


7.8 ±18 


±13 


2.5 < Pt{D*) < 4.5 GeV 


ri{D*) range 


d'^a/drjdpT stat. 


sys. 




[nb/GeV] [%] 


[%] 


-1.5 -1.0 


2.2 ±6 


±9 


-1.0 -0.5 


3.0 ±4 


±9 


-0.5 0.0 


3.6 ±5 


±9 


0.0 0.5 


3.0 ±5 


±9 


0.5 1.0 


2.3 ±7 


±9 


1.0 1.5 


1.8 ±14 


±9 


4.5 < Pt{D*) < 12.5 GeV 


r]{D*) range 


d'^a/drjdpT stat. 


sys. 




[nb/GeV] [%] 


[%] 


-1.5 -1.0 


0.070 ±10 


±12 


-1.0 -0.5 


0.14 ±6 


±11 


-0.5 0.0 


0.22 ±6 


±11 


0.0 0.5 


0.24 ±5 


±11 


0.5 1.0 


0.18 ±6 


±11 


1.0 1.5 


0.11 ±10 


±12 



Table 6: Bin averaged double differential cross sections for inclusive D* production in bins of 
r]{D*) for three ranges in pt{D*) with their statistical and uncorrelated systematic uncertainties. 
The normalisation uncertainty of 6.0% is not included. 
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HI D*-tagged dijet cross sections 


ri{D* jet) range 


da/dri{D* jet) 


Stat. 


sys. 




[nb] 


[%] 


[%] 


-1.5 -1.0 


2.3 


±12 


±11 


-1.0 -0.5 


3.2 


±7 


±8 


-0.5 0.0 


3.9 


±7 


±8 


0.0 0.5 


3.9 


±8 


±8 


0.5 1.0 


3.4 


±9 


±8 


1.0 1.5 


2.8 


±14 


±8 


r] (other jet) range 


dd/dr] (other jet) 


Stat. 


sys. 




[nb] 


[%] 


[%] 


-1.5 -1.0 


1.2 


±15 


±11 


-1.0 -0.5 


1.3 


±13 


±9 


-0.5 0.0 


2.1 


±10 


±8 


0.0 0.5 


2.6 


±9 


±8 


0.5 1.0 


2.7 


±8 


±8 


1.0 1.5 


2.9 


±8 


±8 


1.5 2.2 


2.5 


±10 


±8 


2.2 2.9 


2.2 


±15 


±8 


Pt(D* jet) range 


da/dpT{D*iet) 


Stat. 


sys. 


[GeV] 


[nb/GeV] 


[%] 


[%] 


3.5 5.0 


2.7 


±8 


±8 


5.0 8.0 


1.4 


±5 


±7 


8.0 15.0 


0.17 


±7 


±7 


Pt (other jet) range 


dcr/dpr (other jet) 


Stat. 


sys. 


[GeV] 


[nb/GeV] 


[%] 


[%] 


3.5 5.0 


3.0 


±7 


±8 


5.0 8.0 


1.2 


±5 


±7 


8.0 15.0 


0.24 


±7 


±10 



Table 7: Bin averaged single differential cross sections for D*-tagged dijet production in bins 
of 1] and pt of the D* jet and the other jet with their statistical and uncorrelated systematic 
uncertainties. The normalisation uncertainty of 6.0% is not included. 
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HI I>*-tagged dijet cross sections 


At] range 


da/dArj 


Stat. 


sys. 




[nb] 


[%] 


[%] 


-3.0 -2.0 


0.24 


±33 


±13 


-2.0 -1.0 


0.85 


±12 


±9 


-1.0 0.0 


1.7 


±9 


±8 


0.0 1.0 


2.4 


±7 


±8 


1.0 2.0 


2.5 


±7 


±8 


2.0 3.0 


1.6 


±11 


±8 


3.0 4.0 


0.63 


±21 


±12 


4.0 4.4 


0.22 


±79 


±31 


A(^ range 


da/d\Aip 


Stat. 


sys. 


[deg.] 


[nb/deg.] 


[%] 


[%] 


110 


0.0066 


±24 


±8 


110 150 


0.057 


±8 


±8 


150 170 


0.20 


±5 


±7 


170 180 


0.28 


±6 


±8 


Mx range 


da/dMx 


Stat. 


sys. 


[GeV] 


[nb/GeV] 


[%] 


[%] 


30 75 


0.075 


±4 


±7 


75 120 


0.069 


±7 


±7 


120 250 


0.024 


±11 


±7 


x^ range 


da/dxj 


Stat. 


sys. 




[nb] 


[%] 


[%] 


0.00 0.45 


4.9 


±15 


±9 


0.45 0.75 


11 


±7 


±8 


0.75 1.00 


17 


±4 


±7 



Table 8: Bin averaged single differential cross sections for D*-tagged dijet production in bins 
of Ai], \Aip\, x^ and Mx with their statistical and uncorrelated systematic uncertainties. The 
normalisation uncertainty of 6.0% is not included. 
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HI D*-tagged dijet cross sections 




x^ < 0.75 






Aip range 


da/d|Av9| 


Stat. 


sys. 


[deg.] 


[nb/deg.] 


[%] 


[%] 


110 


0.0057 


±28 


±9 


110 150 


0.040 


±12 


±9 


150 170 


0.10 


±10 


±9 


170 180 


0.12 


±13 


±10 




x^ > 0.75 






A(/9 range 


da/d|Av9| 


Stat. 


sys. 


[deg.] 


[nb/deg.] 


[%] 


[%] 


110 


0.0009 


±34 


±12 


110 150 


0.017 


±11 


±8 


150 170 


0.097 


±6 


±8 


170 180 


0.16 


±6 


±9 



Table 9: Bin averaged single differential cross sections for D*-tagged dijet production in bins 
of \Alp\ in two regions of Xy with their statistical and uncorrected systematic uncertainties. The 
normalisation uncertainty of 6.0% is not included. 
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Figure 1: Distribution of AM for D* candidates a) in the inclusive D* sample and b) in the D* 
tagged dijet sample. The fit function is also shown. 
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Figure 2: Single differential D* cross section as a function of pt{D*), r]{D*), W^p, and z{D*) 
compared to Pythia and Cascade predictions. Here and in the following figures the inner 
error bar depicts the statistical error and the outer shows the statistical, and uncorrected sys- 
tematic and normalisation uncertainty added in quadrature. The normalised ratio R^°^'^ (see 
text) is also shown. 
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Figure 3: Single differential D* cross section as a function of pt{D*), r](D*), Wyp, and z{D*) 
compared to the next-to-leading order predictions of FMNR, GMVFNS and MC@NLO. The 
normalised ratio R^°^^ (see text) is also shown. 
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Figure 4: Double differential D* cross section as a function of r](D*) for three bins oi pt{D* 
compared to PYTHIA and Cascade predictions. 
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Figure 5: Double differential D* cross section as a function of r](D*) for three bins oi pt{D* 
compared to the next-to-leading order predictions of FMNR, GMVFNS and MC@NLO. 
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Figure 6: Single differential cross section for D*-tagged dijet production as a function of pt 
and T] of the D* jet and the other jet compared to Pythia and Cascade predictions. The 
normalised ratio i?"™™ (^^^^ ^^^^^^ j^ ^^^ shown. 
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Figure 7: Single differential cross section for D*-tagged dijet production as a function of pt 
and r] of the D* jet and the other jet compared to MC@NLO predictions. The normalised ratio 

^norm ^^gg j-g-j^j-^ j^ ^^SQ shoWn. 
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Figure 8: Single differential cross section for D*-tagged dijet production as a function of the 
difference in pseudorapidity A?] and in azimuthal angle Ay^ between the other jet and the -D* jet, 
the mass Mx and x^ compared to Pythia and Cascade predictions. The normalised ratio 
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Figure 9: Single differential cross section for D*-tagged dijet production as a function of the 
difference in pseudorapidity A?] and in azimuthal angle A(/9 between the other jet and the -D* jet, 
the mass Mx and x^ compared to MC@NLO predictions. The normalised ratio _R°™™ (see text) 
is also shown. 
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Figure 10: Single differential cross section for D*-tagged dijet production as a function of the 
difference in azimuthal angle A(/3 in two regions of x^ compared to predictions of Pythia, 
Cascade and MC@NLO. The normalised ratio i^-^Q™ (see text) is also shown. 
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